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Preface 

This  report  is  concerned  with  the  effects  of 
vibration  on  convective  heat  transfer  rates.  I  hope  it 
will  add  some  small  contribution  to  the  world  of  ‘know¬ 
ledge  and  be  of  some  use  in  future  studies  of  the  subject. 
The  research  itself  has  been  time  consuming  -yet 
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interesting,  at  times  exasperating  yet  not  without  per¬ 
sonal  satisfaction.  My  one  regret  is  the  lack  of  time 
with  which  to  follow  up  on  the  unexpected. 

I  wish  to  express  my  appreciation  to  Dr.  Andrew  J. 
Shine,  Head,  Department  of  Mechanical  Engineering,  Air 
Force  Institute  of  Technology.  Not  only  did  he  assist  me 
in  selecting  this  particular  project,  but  as  my  thesis 
advisor  he  gave  me  timely  and  invaluable  suggestions. 

Mr.  John  Flahive  and  Mr.  Richard  Brown  provided  vital 
assistance  and  solved  many  pratical  problems  associated 
with  the  assembly  of  equipment.  Finally,  I  wish  to  thank 
Mr.  Millard  Wolfe,  Foreman  of  the  school  shops,  for  his 
assistance  and  interest  without  which  this  project  could 
not  have  been  completed. 


Leon  H.  Chaffee 
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Abstract 

This (  study,  i-s  a  follow-on  investigation  of  the  effects 
of  sinusoidal  vibration  on  the  heat  transfer  rate  from 
cylinders  in  free  convection  in  air.  The  purpose  was  to 
obtain  and  correlate  heat  transfer  data  in  the  region  of 
the  critical  Reynolds  number,  to  investigate  the  effect  of 
the  vibration  at  higher  Reynolds  numbers  through  the  use  of 
a  larger  diameter  cylinder,  and  to  obtain  Schlieren  photo¬ 
graphs  of  the  boundary  layer  in  the  vicinity  of  the  criti¬ 
cal  Reynolds  number.  Three  stainless  steel  cylinders  with 
0.12  in,  0.25  in,  and  0.75  in  diameters  were  vibrated  at 
frequencies  from  10  to  50  cycles  per  second,  at  amplitudes 
from  0.0185  to  0.765  inches,  and  at  surface- to-ambient 
temperature  potentials  of  50F,  100F,  and  200F.  Results 
show  that  each  cylinder  displays  a  similar  characteristic 
pattern  progressing  from  a  region  in  which  the  heat  trans¬ 
fer  rate  is  independent  of  vibration,  through  a  transition 
region,  to  a  region  where  the  heat  transfer  rate  generally 
parallels  the  recommended  forced  convection  curve  of 

McAdams  and  is  a  function  only  of  the  vibration  intensity. 

*>• 

An  increased  frequency  shifts  the  transition  region  in  the 
direction  of  higher  vibration  intensities.  Schlieren 
photographs  show  a  considerable  increase  in  turbulence 
in  the  boundary  layer  through  the  transition  region. 
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I.  Introduction 


Background 

In  a  free  convective  environment,  the  heat  transfer 
rate  from  a  body  to  its  surrounding  is  determined  by  the 
mechanism  in  the  boundary  layer.  In  order  to  increase  the 
heat  transfer  rate,  the  boundary  layer  must  be  altered  by 
reducing  its  thickness  or  increasing  the  transverse  fluid 
movement  in  the  boundary  layer  or  both.  Vibration  is  one 
method  of  producing  increased  transverse  fluid  motion. 

A  knowledge  of  the  influence  of  vibrations  on  the  rate 
of  heat  transfer  by  free  convection  from  a  heated  surface 
would  be  of  considerable  interest  to  the  design  engineer. 
Unfortunately,  to  date  there  is  a  limited  knowledge 
concerning  the  physical  mechanisms  of  the  processes 
involved. 

A  number  of  studies  have  been  conducted  over 
the,  past  two  decades  in  order  to  ascertain  transverse 
vibration  effects  on  the  convective  heat  transfer  rate 
from  both  flat  plates  and  cylinders.  Two  of  the  more 
recent  studies  involved  the  use  of  sinusoidal  vibration, 
a  method  of  vibration  not  used  previously. 
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In  1962,  Shine  and  Jarvis  investigated  the 
sinusoidal  vibration  effects  on  the  convective  heat 
transfer  rate  from  heated  cylinders  in  air.  They 
mounted  the  test  cylinder  horizontally  with  the  ends 
rigidly  clamped.  The  cylinder,  excited  near  one  end, 
was  vibrated  at  its  resonant  frequency  in  a  sinusoidal 
wave  form  in  the  vertical  plane.  Cylinders  of  0.032  in 
and  0.072  in  diameter  were  vibrated  from  15  to  75  cycles 
per  second  and  at  amplitudes  of  from  0.002  to  0.99  in 
(Ref  4:2). 

In  1965,  Watson  conducted  an  extension  of  the 
work  of  Shine  and  Jarvis  to  investigate  the  effects 
of  sinusoidal  vibration  at  higher  vibration  Reynolds 
numbers  (based  on  the  average  vibration  velocity  4af 
and  cylinder  diameter),  at  two  selected  surf ace- to - 
ambient  temperature  potentials,  and  at  different 
positions  along  the  cylinder.  Cylinders  of  0.072  in, 
0.12  in,  and  0.25  in  diameter  were  vibrated  from  16 
to  80  cycles  per  second  and  at  amplitudes  of  from 
0.05  to  1.5  in.  He  obtained  a  maximum  vibration 
Reynolds  number  of  1100  with  the  0.25  in  diameter 
cylinder  (Ref  5:3). 
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From  the  studies  conducted  to  date  the  following 
summary  of  conclusions  were  made  (Ref  3:3,  4:2,  5:2): 

a. .  The  heat  transfer  rate  was  independent 
of  the  direction  of  transverse  vibration  (horizontal 
or  vertical). 

b.  Below  a  critical  vibration  intensity 
(af),  the  heat  transfer  rate  was  unaffected  by  vibration. 

c.  Above  a  critical  vibration  intensity 
the  heat  transfer  rate  increased  with  the  increased 
vibration  intensity  and  was  a  function  only  of  the 
vibration  intensity. 

d.  At  high  vibration  intensities,  the 
effect  of  vibration  was  independent  of  temperature 
and  the  heat  transfer  rate  generally  paralleled 
the  recommended  forced  convection  curve  of  McAdams. 

e.  The  heat  transfer  rate  was  independent 
of  the  position  of  the  test  point  along  the  sinusoidally 
vibrated  cylinder  and  only  dependent  on  the  vibration 
intensity  at  the  position. 

Purpose 

The  purpose  of  this  study  was  to  extend  the 
investigation  of  the  effect  of  sinusoidal  vibration 
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on  the  convective  heat  transfer  rate  from  heated  cylinders 
to  air.  The  specific  objectives  were  (1)  to  obtain  and 
correlate  heat  transfer  data  in  the  region  of  the  critical 
vibration  Reynolds  number  (dimensionless  vibration 
intensity)  with  the  same  size  cylinders  as  Watson  used 5 
(2)  to  investigate  the  effect  at.  higher  Reynolds  numbers 
through  the  use  of  a  larger  diameter  cylinder  and  a 
redesign  of  the  vibrating  apparatus  to  accommodate  the 
larger  cylinder;  and  (3)  to  obtain  Schlieren  photographs 
of  the  boundary  layer  and  surrounding  fluid  in  the  region 
of  the  critical  Reynolds  number. 
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II .  Apparatus 

The  apparatus  used  in  this  study  consisted  basically 
of  (X)  three  electrically  heated  circular  test  cylinders, 

(2)  a  vibrating  assembly  into  which  a  cylinder  was  clamped, 

(3)  a  protective  enclosure  to  seclude,  the  cylinder  from 
random  convective  air  currents,  (4)  a  photographic 
arrangement  used  to  determine  vibration  amplitude,  and 
(5)  a  Schlieren  apparatus. 

Test  Cylinders 

The  test  cylinders  were  42  to  46  in  long  sections  of 
polished,  stainless  steel,  thin-wailed  tubing  with  0.12  in, 

i 

0.25  in,  and  0.75  in  outside  diameters.  Heat  was  provided 
by  passing  adjustable  direct  electrical  current  through 
the  cylinder.  An  internally  attached  iron-constantan 
thermocouple  permitted  instantaneous  monitoring  of  the 
cylinder  wall  temperature  at  the  test  position,  the 
point  of  maximum  vibration  amplitude. 

Vibrating  Assembly 

Two  different  vibrating  assemblies  were  used  during 

% 

the  course  of  this  study.  The  first,  used  with  the  0.12  in 
and  0.25  in  diameter  cylinders,  was  the  same  assembly  used 
by  Watson  (Ref  5:5).  The  cylinder  was  rigidly  clamped  at 
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both  ends.  An  electromagnetic  vibrator  was  positioned 
directly  below  the  cylinder  and  six  inches  from  the  left 
end  clamp.  Variable  frequency  and  amplitude  vibrations 
were  transmitted  to  the  cylinder  through  the  vibrator 
drive  rod  which  also  served  to  constrain  vibratory  motion 
to  the  vertical  plane.  A  schematic  diagram  of  this 
assembly  is  shown  in  Figure  1.  An  adjustable  but  constant 
axial  tension  was  applied  to  the  cylinder  at  its  right  end 
to  permit  selected  variation  in  the  natural  response 
frequency  of  the  cylinder  and  to  compensate  for  thermal 
expansion.  The  second  vibrating  assembly,  used  with  the 
0.25  in  and  0.75  in  diameter  cylinders,  was  designed  to 
accommodate  the  greater  forces  required  to  vibrate  the 
larger  diameter  cylinders.  The  clamps  were  pinned  to 
allow  rotation  in  the  vertical  plane.  The  vibrator  drive 
rod  was  connected  to  a  moment  arm  extension  of  the  right 
cylinder  clamp.  Constraint  of  vibratory  motion  to  the 
vertical  plane  was  effected  both  through  the  pinned  clamp 
arrangement  and  through  increased  clamp  size.  A. 
schematic  diagram  of  this  assembly  is  shown  in  Figure  2. 
Axial  tension  for  this  assembly  was  applied  to  the  left 
end  of  the  cylinder. 
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Enclosure 

The  test  cylinder  with  its  associated  clamps  and 
clamp  supports  of  the  Watson  assembly  was  mounted  in 
a  48  x  12  x  30  in  enclosure.  For  the  improved  assembly, 
the  enclosure  dimension  were  60  x  20  x  30  in.  In  each 
case  the  top  and  bottom  were  vented  to  permit  free 
convection  around  the  cylinder.  An  iron-constantan 
thermocouple  was  located  within  each  enclosure  at 
cylinder  level  to  sense  ambient  temperature. 

Photograph  Arrangement 

A  parallel  light  beam  was  passed  through  a  glass 
windoiv  in  the  front  panel  of  the  enclosure,  across  the  test 
cyx-.iuer,  and  through  a  l/32  in  wide  slit  in  the  rear 
panel.  The  resulting  shadow  of  the  vibrating  test  cylinder 
was  then  directed  onto  a  rotating  mirror  which  swept  it 
across  the  open  aperture  to  a  Polaroid  camera.  The 
resulting  photograph  provided  a  measureable  amplitude  of 
vibration.  A  schematic  diagram  of  the  photographic 
arrangement  is  shown  in  Figure  3.  For  the  Schlieren 
photographs  the  slit  was  removed  and  a  spark  lamp  provided 
the  open  aperture  time  for  the  photograph. 
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Parabolic  Mirror 


Figure  3 

Schematic  of  Photographic  Arrangement 
for  Measuring  Vibration  Amplitude 
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III.  Experimental  Procedure 

The  test  cylinder  was  initially  heated  without 
vibration  to  one  of  the  following  arbitrarily  chosen 
temperature  potentials:  SOP,  100P,  200F.  The  total 
electrical  power  El  delivered  to  the  cylinder  was  then 
determined  by  recording  the  metered  values  of  the  current 
I  and  the  voltage  drop  E.  A  photograph  of  the  static 
cylinder  was  taken  to  determine  the  scale  factor  for 
subsequent  vibrational  amplitude  measurements.  The 
temperature  potential  and  power  measurements  provided 
data  for  determining  the  free  convective  heat  transfer 
rate. 

The  cylinder  was  then  vibrated  at  its  response 
frequency  at  incremental  increases  in  amplitude.  The 
frequency  and  amplitude  measurements  provided  data  for 
determining  the  vibration  Reynolds  number.  This  procedure  > 
was  continued  until  a  decrease  in  temperature  potential  was 
noted  indicating  an  increase  in  the  heat  transfer  rate. 

Prom  this  point  on  the  power  delivered  to  the  cylinder 
was  increased  incrementally  and  the  amplitude  of  vibration 
adjusted  until  the  temperature  potential  had  stablized  at 
the  originally  chosen  value.  Recordings  of  the  power 
delivered  to  the  cylinder,  temperature  potential,  and 


10 


GAM/ME/66B-I0 


frequency  of  vibration  were  made  and  a  photograph  of  the 
vibrating  cylinder  was  taken  for  each  adjusted  change  in 
amplitude.  This  procedure  was  continued  for  a  series  of 
test  runs  during  which  the  temperature  potential  was 
maintained  at  the  originally  chosen  value.  The  double 
amplitude  of  vibration  2a  was  determined  for  each 
photograph  using  a  microscopic  comparator  graduated  to 
0.0001  in.  The  experimental  ranges  were  as  follows: 
frequency  f,  10  to  50  cps;  amplitude  a,  0.0185  to  0.765  in. 

The  maximum  error  in  the  determination  of  the 
Reynolds  number  arising  from  measurements  of  a  and  f  was 
quite  large  at  extremely  low  Reynolds  numbers;  however, 
above  a  Reynolds  number  of  75  the  range  of  error  was 
generally  from  two  to  five  per  cent.  The  maximum  error  in 
the  determination  of  the  heat  transfer  coefficient, 
expressed  as  the  Nusselt  number  Nu,  arising  from  measure¬ 
ments  of  temperature  potential  tw-ta,  I,  and  E  was 
calculated  to  be  7.1  per  cent.  As  the  values  of  Re  and  Nu 
increased  the  associated  error  decreased. 

Schlieren  photographs  of  the  0.75  in  diameter 
cylinder  heated  to  a  temperature  potential  of  200F  were 
taken  at  selected  intervals  through  the  entire  sequence 
of  data  obtained  for  that  temperature  potential. 
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IV.  Results 

The  results  of  this  study  are  presented  in  Figures 

4  through  7.  Heat  tranfer  coefficients  are  expressed  as 
Nusselt  numbers  and  vibration  intensities- are  expressed  as 
Reynolds  numbers. 

There  exist  for  each  cylinder  (1).  a  region  at  low 
vibration  intensities  in  which  vibration  has  no  effect  on 
the  heat  transfer  rate;  (2)  a  region  at  higher  vibration 
intensities  in  which  the  variation  in  the  heat  transfer 
coefficient  generally  parallels  the  recommended  forced 
convection  curve  of  McAdams;  and  (3)  a  characteristic 
transition  region  in  between. 

There  is  conclusive  evidence  that  the  rate  of  heat 
transfer  in  the  transition  region  is  not  a  function  of 
vibration  intensity  alone.  Data  presented  in  Figures 4  and 

5  for  temperature  potentials  of  92F  and  100F  show  that  an 
increase  in  vibration  frequency  shifts  the  transition 
region  in  the  direction  of  higher  vibration  intensities. 

Although  the  variation  of  the  heat  transfer  co¬ 
efficient  with  vibration  intensity  results  in  the  same 
basic  pattern  for  each  cylinder  size  tested,  in  the  region 
where  the  variation  generally  parallels  the  curve  of 
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McAdams  there  does  exist  an  increase  in  displacement  from 
this  curve  with  increased  cylinder  size.  This  is  shown  in 
Figures  4  through  7. 

At  the  higher  vibration  intensities,  the  variation 
in  the  heat  transfer  coefficient  with  vibration  intensity 
for  a  given  cylinder  size  is  independent  of  the  temperature 
potential  and  solely  a  function  of  the  vibration  intensity. 

This  is  shown  in  Figure  6. 

The  Schlieren  photographs  of  the  0.75  in  diameter 
cylinder  at  a  temperature  potential  of  100F  are  presented 
in  Figures  13  through  19. 


GAM/ME/66B-10 


fl±ri  uflArA 

Taillktiu: 


if-**'*'  '  Ht  ~*~*~*  ! 


iflfliflfAA{iA: 


till 

4+  *~t' 

m 

flfl 

flfl 

m 

iffl 

4tii 

Ha 

•flr  H- 

i‘4  i4 

ait 

4fl£ 

flfl 

-Hfl 

Wf 

rttr 

iffl 

(Jit 

■tfl* 

tef 

■M;* 

4{A- 

•♦-tt  t 

flfl 

,4rt 

f  4  +  4 
,4+i*. 

►  f  •  fr 

44  ifl 

4  4-4  4- 

fl?|: 

flfl 

fl 

rr-t-*  »• 

Ait 

ixii 

4  i  ►*. 

rHl 

fl4r 

+•*+* 

*-♦4*. 

1  *-i  4* 
•4*4  4‘ 4 ■ 
•4  *•“  *' 

•  H  f 

*4tf 

.m:| flAJurAiflA  flfl  flA  AfliflA  flfl  flfl  :iu  flit 
flirt- it  tflflttfl  ufliflflflflnAA  flfl  flfl  fltflflA 

■**  •  »-*-  ■•»«.». 4-  -fifl  ■«  4~*  *•*'  *»-»-♦-•-  -4"44-»4  »•«■»>  *-*~**  ■*-*  ♦  » 

flfl.  flfljrrrr:  flfl:  AA  flfl  mflfljrptiiflfli.  i  r:  r  it:  t 

flfl  flit  flfl:  flxr  iflfltfli:  AulAA  fliflflfl  flFPflU. 
ufl flH fljflrfli  flit  jfli iflflflii  flnjflfl  fljflgfli 

-*-4»^  fl  -♦-•-ffl-  -4.»A-4-  ~»A-fA-  •*--»'«-»-  tt  i  t  ^ K* 4~T"*  *  *  t * ** 

Ifl  Hr  fl~  HH  McAdams  ~^Ei 

rRri+H'  fl— r*f  —  -f  44*-  -*•*-»  • «  * .  *W*  y~ 

ti‘***4'  |  •***-*•  -(*■»*•  *  *  * » *  *  |*  ■ «  * 

-fU|.  .*  fAA-  4  *  »  4*  *4  4.  **  j -»  *  *-'+4*«  «  »♦  I  •  %  *  *  0  '  *  * 

►iit tlx * i« »i •< •  •  I *-a *.•*•«  .. .*  **» 

rrrrti^i  mtutr  rr  n  ?  v.  ?  *  r  -  ‘  F  ■ «  nil.il— 

-*-■»- ►  *  i i  »!■«••  >  »•  »  1 1 *-•  *«  J  »«  »  -  »  •  *  *  *  • 

:u:« uu.  flflfl*:: :  urAArpud.:::  ,:*; 


Miaipi 

'jj'j'j'  ^  •j*-r+-  U-||  jfBjjffi 

ilflflfl^Suflj' 

jifllfl  j  rgl  flkflf  ffljfltfl: 

flfl:  flfl  Sfs?  flit  fljflflfl 

flu  Jflgjfl::  flfl.  flfljflH: 


pfl) 


”fl  jufl AflAA rnr  :::: i:::|::;:  uu]:;:: 

1.0  1.4  1.8 

Log  Re 


X?  flflf-i  i-i.  I444  flfl-fH+m- 
flflfflr:  flfl  flflflflflpxrt 

-4  1  ■*  t  4<  *•■  if  1 44  i  M  M4444-}->-»f  f 

fliSflMil 

h  A  ifl ::  IflilflA  iHiHilfl 

»-«♦♦.  «-*»«.  .  -  *  «  i-f  *  i  »  l  4-»flt  I  *  vi-*' 
t  *•  rt  *~*  'f  7 1  f  ^  *  j''' 1 1*;  ;  * 

4.^  *  1*  •  ♦  r  J  i  *  t  i-{  *  i  i  t:l  V4*fl|  {  +444* 
tw-ta  cps  ^jfl 
t  ■  ■ •  4h+ 

•  100F  24-283g 
o  200F  24-28;flir 
■  100F  32-35jtrf 
4  10.QF  4Q~42oxH; 

i  i.flfl  Ml  *  ♦♦HiT+.J  fif 

r:  54} ; ; :  :|:fl  :fl :  rr  j  Iflrjflfl 


Figure  4 

Variation  of  the  Heat  Transfer  Coefficient  with  Vibration 
Intensity  for  the  0.12  in  Diameter  Test  Cylinder 
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Variation  of  the  Heat  Transfer  Coefficient 
with  Vibration  Intensity  at  Selected  Vibration 
Frequencies  for  the  0.25  in  Diameter  Test  Cylinder 
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V.  Discussion  of  Results 

Figure  8  shows  the  variation  of  the  heat  transfer 
coefficient  with  vibration  intensity  at  similar  vibration 
frequencies  for  the  three  cylinders  tested  during  this 
study.  Each  curve  has  a  similar  characteristic  pattern. 
Below  a  critical  vibration  Reynolds  number  the  heat 
transfer  rate  is  independent  of  the  vibration  intensity. 
Above  this  Reynolds  number  the  heat  transfer  rate  increases 
markedly  through  the  transition  region  as  the  vibration 
intensity  is  increased.  At  still  higher  vibration 
intensities  the  heat  transfer  rate  follows  a  curve  that 
generally  parallels  the  recommended  forced  convection 
curve  of  McAdams.  The  displacement  of  the  heat  transfer 
rate  curve  f,rom  the  curve  of  McAdams  in  this  parallel 
region  increases  as  the  cylinder  diameter  is  increased. 

Results  obtained  with  the  0.12  in  diameter  cylinder 
are  compared  in  Figure  9  to  those  obtained  by  Watson  for 
the  same  size  cylinder.  Correlation  is  not  satisfactory 
since  Watson  does  not  show  the  same  characteristic 
transition  pattern.  Examination  of  the  experimental  data 
by  Watson  revealed  that  his  initial  Nusselt  number  (no 
vibration)  is  25  per  cent  higher  than  the  free  convection 
value  recommended  by  McAdams  for  the  cylinder  diameter 
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and  temperature  potential  involved  (Ref  2:176). 
Consequently,  the  data  obtained  by  Watson  is  question¬ 
able.  The  initial  Nusselt  number  obtained  in  this  study 
deviates  from  that  recommended  by  McAdams  by  five  per 
cent. 

In  Figure  10  the  results  of  this  study  with  those 
of  Watson  and  Neely  for  the  0.25  in  diameter  cylinder 
are  compared.  Again,  correlation  with  Watson  is 
unsatisfactory  at  the  lower  vibration  Reynolds  numbers. 

For  this  cylinder,  the  initial  Nusselt  number  he  obtained 
is  19  per  cent  higher  than  the  recommended  McAdams  value; 
therefore,  his  data  is  again  questionable.  The  initial 
Nusselt  number  obtained  in  this  study  deviates  from  the 
recommended  McAdams  value  by  one  per  cent.  Correlation 
with  Neely,  on  the  other  hand,  is  excellent;  the  difference 
in  the  paths  through  the  transition  region  can  be  attrib¬ 
uted  to  different  frequencies  of  vibration.  This  is 
brought  out  more  clearly  in  Figure  11  which  shows  that 
the  transition  region  moves  toward  .lower  vibration 
intensities  as  the  frequency  is  decreased. 

An  additional  significance  of  the  excellent 
correlation  with  Neely  is  the  indication  that  the 
variation  of  the  local  heat  transfer  coefficient  with 
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Comparison  of  Results  for  the 
0.12  in  Diameter  Test  Cylinder 
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vibration  intensity  does  not  depend  upon  the  method  of 
vibration  used.  The  method  of  vibration  used  by  Neely 
was  not  sinusoidal  as  was  used  ’in  this  study.  He  vibrated 
the  entire  cylinder  transversely  in  the  horizontal  plane. 

Figure  12  shows  a  comparison  of  the  test  results 
obtained  with  the  0.75  in  diameter  cylinder  in  this  study 
and  the  0.875  in  diameter  cylinder  by  Fand  and  Kaye 
(Ref  1:495).  Although  direct  correlation  is  not  possible 
since  the  cylinders  are  of  different  size,  the  general 
transition  patterns  appear  similar.  Unfortunately,  Fand 
and  Kaye  were  not  able  to  obtain  data  at  higher  vibration 
intensities  to  show  the  region  generally  paralleling 
the  recommended  curve  of  McAdams.  Fand  and  Kaye  show  a 
slightly  different  curvature  through  the  transition  region 
which  might  be  explained  by  the  fact  that  their  data  was 
obtained  using  higher  and  variable  frequencies,  54  to  225 
cycles  per  second,  rather  than  using  constant  frequencies 
as  in  this  study  (Ref  1:494). 

The  results  of  this  study  show  that  for  a  given 
vibration  frequency  there  is  a  given  path  for  the 
variation  of  the  heat  transfer  coefficient  through  the 
transition  region.  The  supposition  here  is  that  the  over¬ 
all  effect  is  a  function  of  vibration  frequency  when  in 
reality  it  may  be  a  function  of  vibration  amplitude.  In 
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either  case  it  is  apparent  that  the  variation  of  the  heat 
transfer  coefficient  through  the  transition  region  is  not 
soley  a  function  of  the  vibration  intensity,  but  is  indeed 
more  complex.  Additional  qualitative  study  of  this 
phenomenon  is  needed. 

Figure  8  shows  an  increase  in  displacement  from  the 
recommended  curve  of  McAdams  with  increased  cylinder 
diameter.  A  critical  examination  of  the  results  obtained 
in  this  study  and  those  obtained  by  Neely  indicate  that 
the  variation  of  the  heat  -transfer  coefficient  with  vib¬ 
ration  intensity  above  the  transition  region  may  not 
parallel  the  curve  of  McAdams  but  rather  may  be  constantly 
diverging  from  it.  Further  qualitative  study  is  needed  to 
investigate  the  behavior  of  this  possible  divergence, 
particularly  in  the  gap  between  the  0.25  in  and  0.75  in 
diameter  cylinders . 

The  Schlieren  photographs  show  a  dramatic  increase 
in  vortex  turbulence  in  the  boundary  layer  as  the  vibra¬ 
tion  intensity  is  increased  through  the  transition  region. 
This  increase  in  turbulence  provides,  ample  justification 
for  the  rapid  increase  in  the  heat  transfer  rate  with 
increased  vibration  intensity  through  the  transition 
region. 
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Figure  16 

Schlieren  Photographs  at  Position  2 
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Figure  17 

Schlieren  Photographs  at  Position  3 
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Figure  18 

Schlieren  Photographs  at  Position  4 
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Figure  19 

Schlieren  Photographs  at  Position  5 
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VI.  Conclusions  and  Recommendations 
The  results  of  this  study  of  the  free  convective 
heat  transfer  from  vibrating  cylinders  lead  to  the 
following  conclusions: 

1.  There  exist  for  each  cylinder  (1)  a  region 
at  low  vibration  intensities  in  which  vibration  has  no 
effect  on  the  heat  transfer  rate;  (2)  a  region  at  higher 
vibration  intensities  in  which  the  variation  in  the  heat 
transfer  coefficient  generally  parallels  the  recommended 
forced  convection  curve  of  McAdams;  and  (3)  a  character¬ 
istic  transition  region  in  between. 

2,  In  the  region  at  higher  vibration 
intensities  where  the  heat  transfer  rate  curve  generally 
parallels  the  recommended  forced  convection  curve  of 
McAdams,  there  is  an  increased  displacement  from  the 
McAdams  curve  with  increased  cylinder  size  for  the  range 
of  cylinder  sizes  tested.  For  a  given  cylinder  size, 
however,  the  heat  transfer  rate  curve  is  independent  of 
the  temperature  potential  and  dependent  only  on  the 
vibration  intensity  within  the  range  of  temperature 
potentials  used  in  this  study. 
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3.  in  the  transition  region,  the  heat 
transfer  rate  is  hot  solely  a  function  of  the  vibration 
intensity.  An  increase  in  vibration  frequency  shifts 
the  transition  region  in  the  direction  of  higher 
vibration  intensities. 

4.  There  is  a  considerable  increase  in 
turbulence  in  the  boundary  layer  through  the  transition 
region  which  in  turn  justifies  the  marked  increase  in 
the  heat  transfer  rate  through  this  region. 

5.  The  results  of  this  study  are  in 
excellent  agreement  with  the  results  of  Neely  for  the 
0.25  in  diameter  cylinder  (Ref  3:30). 

The  following  recommendations  are  made  for  future 
studies  which  might  be  undertaken: 

1.  That  a  qualitative  study  and  analysis 
of  the  transition  region  phenomenon  be  conducted. 

2.  That  the  divergence  of  the  variation  of 
the  heat  transfer  coefficient  with  vibration  intensity 
from  the  recommended  forced  convection  curve  of  McAdams  be 
investigated. 
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Apparatus 

Arrangement  of  the  apparatus  used  in  this  study  is 
shown  photographically  in  Figures  1  through  4.  These 
photographs  should  supplement  most  of  the  description 
which  follows: 

Test  Cylinders 

The  test  cylinders  were  lengths  of  stainless  steel 
tubing  of  three  different  diameters.  Specific  data  is 


listed  below. 

Diameter 

Length 

Wall  Thickness 

Surface 

0.12  in 

42.6  in 

0.020  in 

Polished 

0.25  in 

43.375  in 

0.028  in 

Polished 

0.75  in 

46.0  in 

0.020  in 

Polished 

An  iron-constantan  thermocouple  was  attached  to  the 
internal  wall  surface  of  each  cylinder  at  the  test  point, 
the  point  of  maximum  vibration  amplitude,  to  provide 
instantaneous  monitoring  of  the  cylinder  wall  temperature. 
The  cylinder  ends  were  closed  with  press-fit  styrofoam  to 
reduce  heat  losses  by  convection  within  the  cylinder. 

Vibration  Assembly  and  Housing 

For  the  Watson  assembly,  the  test  cylinder  was 
mounted  horizontally  and  clamped  rigidly  at  each  end  to 
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a  vertical  post.  An  electromagnetic  vibrator  was 
positioned  directly  below  the  test  cylinder  and  six  inches 
from  the  left  mounting  clamp.  An  aluminum  vibrator  drive 
rod  transmitted  controllable  frequency  and  amplitude 
outputs  to  the  test  cylinder,  and  to  some  degree  con¬ 
strained  the  cylinder  to  vibrate  in  the  vertical  plane. 

An  adjustable  but  constant  axial  tension  was  applied  to  the 
cylinder  at  its  right  end  to  permit  selected  variation  in 
the  natural  response  frequency  of  the  cylinder  and  to 
compensate  for  thermal  expansion.  The  test  cylinder, 
champs,  and  support  posts  were  mounted  in  a  48  x  12  x  30 
in  enclosure  to  provide  protection  from  the  random 
convective  air  currents  in  the  laboratory.  The  enclosure 
was  vented  top  and  bottom  to  permit  free  convection  around 
the  cylinder.  An  iron-constantan  thermocouple  was  mounted 
at  cylinder  level  and  six  inches  from  the  cylinder  test 
point  within  the  enclosure  to  sense  the  ambient  air 
temperature  (Ref  5:5). 

For  the  improved  apparatus,  the  test  cylinder  was 
mounted  horizontally  in  a  three  inch  long  steel  clamp  at 
either  end.  Each  clamp  was  pinned  to  a  nine  inch  high 
clamp  support  which  allowed  rotation  in  the  vertical  plane 
and  constrained  cylinder  vibration  to  the  same  plane.  The 
clamp  supports  were  bolted  to  a  five  foot  long  piece  of 
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5  x  1-7/8  in  inverted  channel  iron  which  inturn  was  bolted 
to  the  top  of  a  26  x  28  x  49  in  cement  pedestal  in  the 
M.  E,  Laboratory.  The  left  clamp  support  was  adjustable 
along  the  channel  iron  to  accommodate  different  cylinder 
lengths.  An  adjustable  but  constant  axial  tension  was 
applied  to  the  test  cylinder  at  the  left  mounting  clamp. 

The  electromagnetic  vibrator  was  bolted  to  the  right  side 
of  the  cement  pedestal.  The  steel  vibrator  drive  ro  was 
connected  to  a  moment  arm  extension  of  the  right  mounting 
clamp.  The  channel  iron,  clamps  and  supports,  and  test 
cylinder  were  mounted  in  a  60  x  20  x  30  in  enclosure  with 
removable  front  and  back  panels.  An  iron-constantan 
thermocouple  for  measuring  the  ambient  temperature,  was 
mounted  at  cylinder  level  and  nine  inches  from  the  cylinder 
test  point. 

Cylinder  Heating  Power 

A  28  volt,  100  amp  power  source  was  available  from 
the  M.  E.  Laboratory  DC  rectifier.  This  source  was  con¬ 
nected  to  a  high  resistance  load  bank  which  permitted  6  to 
8  amp  incremental  changes  in  current.  A  5.8  ohm,  10  amp, 
and  a  2.3  ohm,  4  amp  slide  resistor  were  connected  in 
parallel  with  the  load  bank  to  provide  fine  power  adjust¬ 
ments.  Power  leads  connected  to  metal  clamps  at  either  end 
of  the  test  cylinder  placed  the  cylinder  in  series  with  the 
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resistance  and  an  ammeter.  The  voltage  pick-off  leads 
were  connected  to  the  power  clamps  for  the  0.12  in  and  0.25 
in  diameter  cylinders  and  were  soldered  to  the  0.75  in 
diameter  cylinder  immediately  inside  the  mounting  clamps. 
The  location  of  these  voltage  leads  determined  the  cylinder 
lengths  L  used  in  computing  the  cylinder  surface  area  Aw. 


These  lengths  are 

listed  below: 

Diameter 

Run  Series 

Length 

0.12  in 

A,  B,  C 

40.0  in 

0.25  in 

D,  E,  F 

41.8  in 

0.25  in 

G,  H 

34.6  in 

0.75  in 

I,  J 

37.125  in 

Temperature  Measurement 

Two  temperature  recorders  were  used  during  this 
study.  A  Minneapolis  Honeywell  Brown  Electronik  Recorder 
with  a  50F  to  300F  range  and  12  channels  was  used  for  the 
A  through  H  series  of  runs.  A  Bristol  Dynamaster  Recorder 
with  a  OF  to  200F  range  and  16  channels  was  used  for  the 
I  and  J  series  of  runs.  For  each  recorder  the  ambient  air 
sensing  thermocouple  was  connected  to  one  channel  to 
provide  a  temperature  monitor  for  each  minute  of  operation. 
The  cylinder  wall  thermocouple  was  connected  in  series  to 
each  of  the  remaining  channels  to  provide  a  temperature 
monitor  every  four  to  five  seconds. 
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Vibration  Control  and  Measurement 

Vibration  frequency  was  controlled  through  adjustment 
of  an  audio  frequency  oscillator  contained  in  the  vibrator 
power  supply  cabinet.  Actual  frequency  measurements, 
however,  were  determined  with  a  General  Radio  Company 
strobotaci 

Vibration  amplitude  was  controlled  by  a  power 
rheostat  in  the  vibrator  power  supply.  In  order  to  measure 
the  amplitude,  a  photographic  arrangement  was  used.  The 
light  from  a  mercury  lamp,  unit  was  passed  through  a 
condensing  lens  onto  a  1/32  in  pin  hole.  The  resulting 
point  source  of  light  was  located  at  the  focal  length  of 
a  7.5  in  parabolic  mirror  which  directed  a  parallel  light 
beam  through  a  vertical  1/32  in  slit  in  the  rear  panel  of 
the  enclosure,  across  the  vibrating  cylinder,  and  through 
a  glass  window  in  the  front  panel.  The  shadow  image  of 
the  vibrating  cylinder  was  reflected  by  a  second  parabolic 
mirror  onto  a  2  in  rotating  flat  mirror  which  swept  the 
image  across  the  aperture  to  a  Polaroid  camera  containing 
type  47  film.  The  camera  was  mounted  on  the  end  of  a 
plywood  housing  with  an  extendable  bellows  section  which 
contained  the  aperture  and  a  manually  operated  shutter. 
Knife  scratches  were  inscribed  on  the  photograph  across 
the  vibration  peaks  and  the  double  amplitude  was  measured 


6 


GAM/ME/66B-10 


with  a  microscopic  comparator.  Figure  5  shows  examples 
of  the  photographs  used  to  determine  vibration  amplitude. 

For  the  Schlieren  photographs,  the  rotating  mirror 
was  replaced  with  a  4  in  plane  mirror  and  the  1/32  in  slit 
was  removed  to  provide  a  more  lateral  view  of  the  cylinder. 
A  knife  edge  was  located  between  the  4  in  plane  mirror  and 
the  camera  aperture.  A  five  microsecond  spark  lamp  was 
installed  in  place  of  the  mercury  lamp;  the  spark  feature 
also  served  as  the  camera  shutter. 
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Figure  5 

Photographs  of  Static  and  Vibrating 
0.75  in  Diameter  Cylinder 
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Experimental  Procedure 

The  static  test  cylinder  was  initially  heated  to  a 
selected  temperature  potential.  The  current,  voltage, 
ambient  air  temperature,  and  cylinder  wall  temperature 
were  recorded.  A  photograph  of  the  static  cylinder  was 
taken  to  provide  a  scale  factor  for  computing  subsequent 
vibration  amplitudes.  The  photographic  procedure  consisted 
of  the  following:  (1)  the  room  lights  were  turned  off;  (2) 
the  camera  slide  was  removed;  (3)  as  the  image  from  the 
rotating  mirror  passed  the  camera  aperture,  the  manually 
operated  shutter  was  opened  momentarily;  (4)  the  camera 
slide  was  re-installed;  and  (5)  the  room  lights  turned 
back  on.  This  procedure  was  also  used  for  all  photographs 
taken  of  the  vibrating  cylinder. 

The  vibrator  was  turned  on  and  the  frequency  adjusted 
to  the  natural  response  frequency  of  the  cylinder.  The 
vibration  amplitude  was  increased  incrementally.  The 
current,  voltage,  ambient  air  temperature,  cylinder  wall 
temperature,  and  frequency  were  recorded  and  a  photograph 
was  taken  for  each  incremental  increase  in  amplitude.  This 
procedure  was  continued  until  a  decrease  in  the  temperature 
potential  which  indicated  an  increase  in  heat  transfer 
rate,  was  noted.  From  that  point  on,  the  current  was 
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increased  incrementally  and  the  vibration  amplitude  was 
adjusted  until  the  selected  temperature  potential  was  again 
obtained  and  stabilized.  The  current,  voltage,  ambient 
air  temperature,  cylinder  wall  temperature,  and  frequency 
were  recorded,  and  a  photograph  taken  for  each  incremental 
increase  in  current.  This  procedure  was  continued  through 
each  series  of  test  runs. 
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Development  of  Equations 
Equation  for  Calculating  Nu 

Development  of  the  working  equation  for  determining 
the  local  Nu  at  the  test  point  on  the  cylinder  involved 
the  following  assumptions: 

a.  The  electrical  power  per  unit  volume 
delivered  to  the  cylinder  was  constant  between  voltage 
pick-off  points, 

b.  Axial  heat  conduction  along  the  cylinder 
at  the  test  point  was  negligble. 

c.  Radiation  heat  losses  were  a  function  of 
the  temperature  potential  and  not  a  function  of  vibration. 

Axial  temperature  profiles  were  experimentally 
determined  for  the  0.25  in  diameter  cylinder  at  tw-ta  " 
100F  and  for  the  0.75  in  diameter  cylinder  at  tw-ta  -  50F 
and  tw-ta  =  100F  by  moving  an  iron-constantan  thermocouple 
at  measured  intervals  along  the  inner  wall  of  the  static 
cylinder.  The  resulting  temperature  profiles  are  shown  in 
Figures  6,  7,  and  8.  Since  the  temperature  gradients 
were  essentially  flat  over  a  considerable  portion  of  the 
cylinder  including  the  center  test  point,  the  axial  heat 
conduction  was  shown  to  be  negligible.  Similar  results 
were  arbitrarily  assumed  for  the  0.12  in  diameter  cylinder. 
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The  power  per  unit  volume  was  experimentally 
determined  over  selected  intervals  of  the  static  0.75  in 
diameter  cylinder  at  t  -t  =  50F.  The  initial  interval 

W  3, 

was  8  in,  4  in  either  side  of  the  center  test  point.  The 
second  interval  was  10  in,  5  in  either  side  of  the  center 
test  point.  Subsequent  intervals  were  determined  in  the 
same  manner  across  the  length  of  the  cylinder  until  the 
vibration  test  voltage  pick-off  points  were  reached.  At 
each  interval  the  voltage  drop  across  the  interval  was 
determined  using  a  Weston  Analyzer  with  a  1.6  volt  scale. 
The  current  was  maintained  constant.  Volumes  were 
computed  assuming  a  constant  cross  sectional  area  and 
using  the  selected  interval.  Results  are  shown  in 
Figure  4  and  within  the  accuracy  of  measurements  the 
power  per  unit  volume  was  constant  over  each  interval 
tested.  Similar  results  were  arbitarily  assumed  for 
the  other  cylinders  and  temperature  potentials. 

The  basic  equation  for  the  local  convective  heat 


loss  can  be  given  by 


Q~dV 

—  =  hdA(tw-ta) 


This  can  be  rewritten  as 


Qcftxg*  =  hcdx(tw-ta) 
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Where  c  =  perimeter  of  the  cylinder.  Furthermore, 

since  Ax  ,  1 
V  "  L 

J?  =  hc(tw-ta) 

Qc  -  hcL(  "tw-ta) 

=  ^w(^w"^a) 

h  =  Qc _ 

^w( ^w"^a) 

Nu  =  .  Qr  D  ■ 

AWkf{tw-ta) 

Now  Qc  =  El  -  Qr  where  El,  the  power  inpu%  can  be 
measured  and  Qr  must  be  determined  independently.  Hence, 
the  local  Nu  can  be  determined  using  measurable  quantities 
with  the  exception  of  Qr. 

The  working  equation  with  Qc  expressed  in  watts  is 

QcD 

Nu  =  3.413  Awkf(tw-ta)  (1) 

Equation  for  Calculation  of  Qr 

The  basic  equation  for  calculation  of  Qr  (B/hr)  is 
given  by 

Qr  =  CT-evV-Ta4) 

Evaluation  of  the  emissivity  was  accomplished  by  forcing 
the  static  free  convective  heat  transfer  rate  obtained  by 
the  apparatus  of  this  study  to  agree  with  the  recommended 
free  convection  curve  of  McAdams  (Ref  2:176).  Knowing  Qc 
in  this  case,  Qr  could  be  determined  by  subtracting  Qc 
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from  the  measured  El;  hence,  the  emissivity  could  be 
determined  and  assuming  it  to  be  constant  in  the  temper¬ 
ature  ranges  tested,  it  could  be  reapplied  to  all  data. 
Below  is  a  summary  of  calculations  made  to  determine  the 
emissivity  for  each  cylinder  based  on  the  recommended  free 
convection  curve  of  McAdams. 


Cylinder 

Run 

GrPr 

Nu 

Qr (watts) 

€ 

0.12  in  OD 

A-l 

105 

2.12 

0 

o 

• 

o 

B-l 

143 

2.16 

2.3 

0.22 

0.25  in  OD 

D-l 

919 

3.07 

1.22 

0.139 

E-l 

1226 

3.22 

3.2 

0.138 

G-l 

824 

2.96 

0.94 

0.138 

0.75  in  OD 

1-1 

15850 

5.63 

0.16 

0.015 

J-l 

25000 

6.03 

1.44 

0.059 

The  following  emissivities  were  selected  based  on  the 
above  data: 


Cylinder 

€ 

0.12  in  OD 

0.10 

0.25  in  OD 

0.14 

0.75  in  OD 

0.05 

The  chosen  emissivities  were  used  to  calculate  the 

i 

test  Nu  in  all  data  taken.  Below  is  summarized  the 
deviation  in  the  free  convective  Nu  as  compared  to  McAdams: 
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Cylinder  Run  Nu(tesl 


Nu (McAdams 


Error  % 


0.12  in  0D  A-l  1.91 


B-l  2.27 


2.12 


2.16 


0.25  in  OD  D-l  3.10 


3.07 


E-l  3.20 


3.22 


G-l  2.96 


2.96 


0.75  in  OD  1-1  5.46 


5.63 


J-l  6.07 


6.03 


The  working  equation  for  determining  Qr  in  watts  is 


Qr  -  IZl5Z 


(0.173)  A„, 


■'e  IM 


where  6  is  one  of  the  selected  values  above  depending  on 
the  cylinder  as  identified  by  the  diameter. 

Equation  for  Calculating  Re 

The  average  vibrational  velocity,  4af,  was  used  in 
calculating  Re.  The  working  equation  is 


Re  =  4afD  »  2afD 

12  Vf  6  Uf 
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Temperature  Profile  in  0.25  in  Diameter  Cylinder  with  (t. 
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Error  Analysis 


Electrical  Measurements 

Measurements  of  the  voltage  E  and,  current  I  were 
used  to  calculate  the  total  power  input  to  the  cylinder. 

The  voltmeters  and  ammeters  were  calibrated  by  the 
Electrical  Engineering  Department  against  a  meter  standard 
accurate  to  ^  of  1%.  Instrument  reading  errors  were 
judged  to  be  +  10%  of  the  smallest  scale  reading.  Accuracy 
data  for  the  meters  used  is  presented  in  Table  C-l. 

Frequency  Measurement 

The  vibration  frequency  was  measured  with  a  stro- 
botac  calibrated  prior  to  each  series  of  runs.  The 
accuracy  of  the  instrument  when  calibrated  was  +  1%  of  the 
high  scale.  This  amounted  to  a  possible  error  of  +  0.617 
cps. 

Amplitude  Measurement 

The  double  amplitude  of  vibration  was  measured  with 
a  micro  comparator  between  knife  scratches  scribed  across 
the  amplitude  peaks  on  the  Polariod  film.  The  comparator 
could  be  read  to  the  fourth  decimal  place.  The  primary 
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Electrical  Measurement  Data  Errors 
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source  of  error  was  the  scribing  of  the  double  amplitude. 
This  error  was  estimated  to  be  +  0.001  in.  In  addition, 
the  scale  factor  had  an  estimated  +  0.001  in  error,  The 
maximum  total  error  in  measuring  the  double  amplitude 
was  +  0.002  in. 

Temperature  Measurement 

The  Brown  Recorder,  used  in  the  A  through  H  series 
of  runs,  was  calibrated  over  its  50F  to  300F  range  by 
personnel  in  the  Aerospace  Research  Laboratories.  The 
entire  temperature  measuring  system,  recorder  and  iron- 
constantan  thermocouples,  was  re-checked  at  ambient  air 
and  boiling  water  temperatures  against  a  mercury  in 
glass  thermometer.  Agreement  was  within  the  ability  to 
read  the  scales.  The  Bristol  Dynamaster  Recorder,  used 
in  the  I  and  J  series  of  runs,  was  calibrated  over  its 
OF  to  200F  range  with  an  accurate  Leeds  and  Northrup 
millivolt  potentiometer.  Since  for  either  recorder  the 
manufacturer  could  not  guarantee  an  accuracy  of  better 
than  +0.03  millivolts  in  the  mechanical  slide  linkage, 
a  possible  error  of  +  IF  was  assumed. 


GAM/ME/66B-10 


Overall  Error 

The  maximum  possible  error  in  Nu  is  given  by 
ANu  - 

The  maximum  possible  error  in  Re  is  given  by 


3 . 413D  I"  E  A  I  IAE  El  A ( t,„ - 1 
Awkf  (tw-ta)  (inv-^a)  (^vv^a) 


(4) 


ARe  =  — ZiL  (a  A  f  +  f  A  a) 
I2  Vf 


(5) 


The  largest  overall  error  in  Nu  occurred  for  Run  A-l 
and  was  computed  to  be  7,1%.  As  the  value  of  the  Nu 
increased  its  relative  error  decreased. 

The  largest  overall  error  in  Re  occurred  for  Run  A-7 
which  also  had  the  lowest  Re  (=8).  It  computed  to 
This  resulted  as  a  consequence  of  the  relative  values 
of  the  small  measured  amplitude  of  vibration  and  the 
estimated  error  involved  in  this  measurement.  The  error 
in  Re,  however,  dimenished  rapidly  as  Re  was  increased 
and  the  general  range  of  the  Re  errors  was  2%  to  5%. 
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Sample  Calculations 

The  following  calculations  are  based  on  Run  1-34 
for  the  0.75  in  diameter  cylinder.  The  cylinder 
dimensions  are: 

D  =  0.0626  ft  L  =  3.09  ft 
The  recorded  test  data  are: 

ta  =  87F  I  -  49.0  amps  2a  -  0.6430  in 

txv  =  137F  E  ■  1.26  volts  f  =  43.2  cps 

The  following  computations  were  performed: 
tf  =  tA  -  137  +  87  H2F 

x  2  2 

tw-ta  =  137-87  =  50F 

El  =  (1.26)(49.0)  =  61.74  watts 
Aw  =  7TDL  =  ( 3 . 1416 ) ( 0 . 0626 ) ( 3 . 09 )  =  0.609  ft2 
The  fluid  properties  at  tf  =  112F  were  determined 
to  be  (Ref  1:504): 

Uf  -  18.83x10-5  ft2/sec 
kf  =  0.01594  B/hr  ft.  F 
Pr  =  0.704 

The  radiation  heat  loss  was  computed  from  Equation 


Qr  =  ^^12(0.173)^6 

Tw 

4 

Taf 

60  w 

100 

i°9l  J 

where  the  emissivity  €  =  0.05  (from  Appendix  B)  and  Qr  is 
given  in  watts . 
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ANu 


3.413D  f  EAI  +  IAE  El  A(tw-t-) 

Awkf  (tw"ta)  (tw_ta)  (tw_ta) 


At  =  0.475  amps  AE  =  0.0095  volts  A(tw-ta)  -  IE 

( 1 .26X0.475 )  (49. CM 0.0095).  (61.74)(1) 
50  +‘  50  2500 

ANu  =  1.01 


ANu  =  (3.413X0.06261 
(0.609)(0. 01594) 


Error  in  Nu  =  3.76% 

The  overall  error  in  the  Reynolds  number  from 

Equation  (5)  is 
4D 

ARe  =  -j2^7(aAf  +  f  Aa) 

Af  =  0.617  cps  Aa  =  0.001  in 

ARe  =  (i2)fiitS^)[(0-321s)(0-617)+  (43-2)(0-001)] 

ARe  =  24.4 


Error  in  Re  =  1.59% 
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